Multichannel cold collisions between metastable Sr atoms 
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We present a multichannel scattering calculation of elastic and inelastic cold collisions between 
two low-field seeking, metastable **Sr [{5s5p)^ P2] atoms in the presence of an external magnetic 
field. The scattering physics is governed by strong anisotropic long-range interactions, which lead 
to pronounced coupling among the partial waves of relative motion. As a result, nonadiabatic 
transitions are shown to trigger a high rate of inelastic losses. At relatively high energies, T > 100/^K, 
the total inelastic collision rate is comparable with the elastic rate. However, at lower coUisional 
energy, the elastic rate decreases, and at T ~ Ia^K, it becomes substantially smaller than the inelastic 
rate. Our study suggests that magnetic trapping and evaporative cooling of **Sr [{5s5ip)'^ P2] atoms, 
as well as ^"Ca [{As'ip)^ P2], in low-field seeking states will prove difficult to achieve experimentally. 

PACS numbers: 34.50.-s,34.20.Cf, 34.20.Mq 



In recent years, the trapping and cooling of alkaline- 
earth atoms has been pursued with increased interest. 
Trapping experiments with ^^Sr 0, S H S and -^^Ca 
[6|,|7(| atoms have been carried out, and experiments are 
underway in a number of laboratories. One major thrust 
of this research track is to extend the small class of 
atomic species that have been Bose-Einstein condensed. 
Also, the most prevalent isotopes have closed-shell nuclei, 
which eliminates complications from hypcrfinc structure 
and simplifies the interpretation of cold collision data Q . 

The usual route to quantum degeneracy utilizes mag- 
netic trapping and evaporative cooling of precooled 
atoms. Since ground-state alkaline-earth atoms possess 
no magnetic dipole moment, some experimental groups 
work with a system of alkaline-earth atoms in their lowest 
{nsnpfP2 metastable state 0,11111111. The radiative 
lifetime of this state is of the order of ten minutes ^3 . 
However, as we show in this Letter, it is not spontaneous 
emission that ultimately limits the magnetic-trap lifetime 
of metastable alkaline-earth atoms. Their nonspherical 
electronic structure gives rise to strong anisotropic in- 
teratomic interactions not present in ground-state alkali 
metals, hydrogen, or metastable helium. In the present 
study, these interactions are predicted to generate severe 
inelastic losses. 

Derevianko et al. |^ were the first to theoretically ap- 
proach the intriguing problem of anisotropic long-range 
interactions between metastable alkaline-earth atoms in 
an external magnetic field. Their study revealed impor- 
tant insights, notably the existence of a long-range molec- 
ular potential well. They also provided a first estimate 
of inelastic losses within a two-channel model, but at a 
qualitative level only. Recently, Ref. ^2 performed a 
detailed tensorial analysis of the anisotropic long-range 
forces between two nonspherical atoms in a magnetic field 
and investigated the associated coupling between elec- 
tronic structure and molecular rotational motion. A two- 
body potential was derived that is ideal for multichannel- 
scattering calculations. The present work has the goal 
to carry the aforementioned theoretical studies to the 



next level, and thus to provide guidance to experimen- 
talists. Specifically, we develop a quantitative quantum 
description of the two-body elastic and inelastic collisions 
between alkaline-earth atoms in the (nsnp)^ P2 excited 
state. We focus on strontium atoms but briefiy address 
implications of the present study to Ca and other similar 
atoms. 

Our treatment is based on the multichannel potential 
derived in Ref. [T^ : 

Hi ^ J2 V2KTl[CK®VTKV]a,,, (2) 

X=0,2,4 

which is a function of internuclear distance R. L is the 
operator of the angular momentum associated with the 
rotation of the molecular frame, miod is the reduced mass 
of the system, /x^^^ and /x^^^ are magnetic dipole opera- 
tors of the two atoms, and B is the vector of the magnetic 
field. The term Hi is the effective interatomic interac- 
tion Hamiltonian acting in the electronic subspace of two 
{nsnpY P2 atoms, as indicated by the projection operator 
V . -ffi is a totally invariant tensor, consisting of three dis- 
tinct contributions, which derive from coupling of rank- 
K tensor operators Ck and Tk- Ck is associated with 
molecular rotational motion, Tk describes specific elec- 
tronic couplings. Electric dipole-dipole dispersion forces 
contribute to K — 0, 2, and 4, magnetic dipole-dipole in- 
teraction enters T2, and electric quadrupole-quadrupole 
interaction contributes to T4 . For the incident chan- 
nels of interest, and for magnetic fields of a few gauss 
or stronger, all collision dynamics at cold temperatures 
takes place at interatomic distances of about 100 Bohr 
radii. Thus, using a pure long-range Hamiltonian is jus- 
tified. In the present calculation the potential V is rep- 
resented by a matrix in the coupled angular momentum 
basis and evaluated using Wigner-Racah algebra. 

The potential matrix of Eq. ^ can be diagonalized at 
different R, giving a set of adiabatic potential curves. 
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At large internuclear distances, the potential curves con- 
verge to nine different dissociation limits, corresponding 
to nine different values {M^i = —4, . . . , 4) of the total 
electronic angular momentum projection on the direction 
of the magnetic field. We assume that the initial relative 
kinetic energy is low enough (< lOOfiK), so that we may 
limit our treatment to the s-wave (L ^ 0) initial chan- 
nel. However, other channels with L > are present in 
our model as possible escape channels. The technique of 
magnetic trapping is limited to low-field seeking atomic 
states, i.e., the states with a negative projection of the 
magnetic atomic moment on the direction of the field. 
We assume that initially, before a collision, the two atoms 
have a specific Mtot projection of the total angular mo- 
mentum including molecular frame rotation. The states 
with different Altot projections are not coupled and can 
be treated separately. Thus, Mtot, unlike M^i, is con- 
served during the collision. In the initial channel L — 0, 
thus M„i = Mtnt for this channel. 
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FIG. 1: Adiabatic potential curves of 812^^2] calculated for 
a magnetic field B = 500 G. The partial 
wave taken into account for these curves is 22 . The incident 
channel [M^i =4, L = 0) is marked by a big arrow. The 
avoided crossing, marked by a circle, produces a large inelastic 
rate. Numbers at the bottom of the plot specify the angular 
momentum L of different partial wave channels converging to 
the Mci = 3 Zeeman level. Note, that at distances of about 
100 a.u. different partial waves are mixed. The inset shows 
the same potential curves on a larger energy scale. In the 
inset one can see all nine dissociative limits corresponding to 
nine Zeeman levels of the two separated strontium atoms. 

Figure ^ gives an example of the adiabatic potential 
curves for Mtot = 4 at the magnetic field B = 500 gauss. 
In the present calculation we have included partial waves 
of rotating Sr2 molecular frame up to L = 22. The en- 
trance s-wave channel is marked by a thick line. This 
channel has a potential well situated around R — 90 
a.u. The well is formed by the interplay of Zeeman, 
quadrupole-quadrupole and dispersion interactions be- 
tween the two atoms [H [13. The shape and depth of 
the well changes with the magnetic field. This fact allows 



to change scattering properties of two-body collisions in 
presence of the field. References 0, suggest that at 
some magnetic field the elastic scattering length can be 
made very large, allowing fast evaporative cooling. 

The key feature in Fig. ^ is a pronounced avoided 
crossing at 80-90 a.u. This avoided crossing is mostly 
formed between the entrance-channel curve and a curve 
converging to a lower dissociation limit with {M^i = 
3, L = 14) (see the encircled region in Fig. Other 
states are mixed as well in this region, but they play a 
secondary role. According to the familiar Landau-Zener 
model 1131, such an avoided crossing produces a large 
transition probability from the entrance channel into the 
{Mtot = 3, L = 14) channel. This is because the en- 
ergy of the avoided crossing is close to the total energy 
of the system, and the two curves are almost parallel near 
the avoided crossing. Thus, at low collision energies, one 
would expect the inelastic cross-sections to be relatively 
high. 

In the previous theoretical study [l3|, Derevianko et 
al. have calculated the elastic collisional rates and esti- 
mated the inelastic rates as well. Their conclusion was 
that inelastic collisions do not significantly contribute to 
the total rate. However, in that study, several simplifying 
assumptions were made. The most important is that the 
authors estimated the inelastic losses using a two-channel 
model. The model includes an s-wave channel coupled 
to a g-wave channel by the quadrupole-quadrupole in- 
teraction. Our calculations show, however, that higher, 
L > 6, partial waves are particularly important in the 
loss process. We have found that L = A does not play 
any important role in the losses, which is in agreement 
with Ref. 2l|- This effect of higher L is apparent in 
Fig. n where the channels (Me/ —4, L = 0) and 
(A/gj =3, L' = 14) are strongly mixed around i? = 85 
a.u. At different magnetic field, different partial waves 
play a dominant role in the inelastic collisions. In other 
words, the avoided crossing in Fig. ^can be formed by 
the states (M^i = 4, L = 0) and (M^, = 3, L') with L' 
depending on the field. The inclusion of the higher partial 
waves not only produces larger non-adiabatic dynamical 
effects, but it also modifies the shape of adiabatic curves. 
This would mean that Derevianko et al. [Tlj might un- 
derestimate the inelastic rates by neglecting the channels 
that 1) modify the entrance adiabatic channel and that 
2) control the inelastic losses associated with the strong 
nonadiabatic coupling. The strong multichannel coupling 
not only causes the large inelastic losses; it significantly 
modifies the elastic scattering. This makes the concept 
of adiabaticity, adopted in Ref. |0 , inapplicable to the 
present case. 

The above analysis of the adiabatic potential curves 
develops a qualitative description of two-body collisions. 
To obtain quantitative results, we solve the coupled ra- 
dial Schrodinger equations for the potential operator of 
Eq. n 
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FIG. 2: Elastic (upper panel) and total inelastic (lower 
panel) rates of the Sr[('^P2)]+Sr[(^P2)] collisions as a func- 
tion of magnetic field for several incident relative kinetic en- 
ergies of the two strontium atoms. The entrance channel is 
{Mel = 4, L = 0). At 1 mK, incoming partial waves with 
L > can penetrate their respective angular momentum bar- 
riers. Their contribution is not included in this plot. (See 
Fig. H 



Here E is the total energy of the system, *t(i?) is the mul- 
ticomponent radial wavefunction, / is the identity opera- 
tor. When solving the equation, we work in the coupled 
channel (diabatic) representation rather than in the adi- 
abatic representation. This choice is justified by the fact 
that in the diabatic representation, all diagonal and non- 
diagonal couplings are smooth functions of R. Out of 
the possible diabatic basis sets we adopt the atomic set, 
which diagonalizes the two-body potential V" at i? = oo. 

For a given asymptotic coUisional energy in the en- 
trance channel, we obtain the total scattering matr ix S 
using the one-dimensional i?-matrix approach 
The i?-matrix is obtained solving the Schrodinger equa- 
tion within R < Rq — 1500 a.u. using finite elements 
[isjl as a representation basis. Having obtained the total 
S^-matrix, the elastic and inelastic rates ai as function of 
E are given by familiar formulas (in atomic units): 



ki' n 

CTi.i' = T 



(4) 



where k'^ — i/2mred(-£' ^ Ei> ) is the incident de Broglie 
wavenumber, Eii is the energy of the incident channel, 
and ai^i' is the corresponding cross-section. The entrance 



channel is specified by index i' . In our discussion of colli- 
sion rates, we do not utilize the scattering length, because 
it is not energy-independent at energies E > lOO/iif for 
this system. 

In order to check the reliability of the present results, 
we tested convergence with respect to different numerical 
and physical parameters of the problem. Two of these, 
in particular, deserve mention here. The first is a conver- 
gence test of the calculated 5-matrix with respect to the 
position of the inner boundary in the numerical calcula- 
tions. Since we have neglected the short-range molecular 
interaction (we do not have the corresponding data), we 
have to be sure that no flux penetrates into the inner 
region, i.e. the coUisional dynamics is controlled only by 
long-range potentials. As a result, we have found that 
for the (Me/ = 4) entrance channel, at low energies be- 
low ImK, no flux penetrates into the short-range region. 
However, the calculation for all the other low-field seek- 
ing magnetic states, {M^i = 1,2,3), shows that there is 
a substantial probability of penetration, and therefore, 
the inelastic rate is expected to be large. These entrance 
channels cannot be treated unless one includes the short- 
range molecular dynamics. Thus, the following treat- 
ment is concentrated on the {M^i =4, L — 0) entrance 
channel solely. The second test is an assessment of con- 
vergence with respect to the maximum number of par- 
tial waves that are included in the calculations. As was 
discussed above, most of the inelastic losses occur due 
to the avoided crossing in the potential well of the inci- 
dent channel. At the largest magnetic field considered 
here, B=1000 G, the avoided crossing is formed by an 
interaction with the M^i —3, L — 14 channel. Perform- 
ing the full S^-matrix calculation with maximum number 
of partial waves, Lmax, up to 18, we have found that, 
indeed, the inelastic and elastic rates are converged at 
Lmax = 14. This test underlines the effect of the molec- 
ular frame rotation on the scattering process. 

Figure El shows the inelastic and elastic rates as a func- 
tion of the external magnetic field for several incident ki- 
netic energies. The resonant-like behavior of the curves 
around 400 G is caused by the emergence of a quasi- 
bound level in the entrance channel. Due to the strong 
inelastic loss, the resonance is significantly washed out. 
Derevianko et al. ^Jl have predicted a resonance at a 
field of about 1000 G. Since they neglected the inelastic 
losses in their model, their elastic rate diverges to infinity 
at the resonant i?-field. In experiments aimed to obtain 
a Bose-Einstein condensate, one is particularly interested 
in the ratio between the elastic and inelastic rates. Fig- 
ure |21 suggests that for strontium in the {5s5p)^P2 state, 
evaporative cooling will not be efficient at temperatures 
below 10 /iK. Figure O shows the rates as a function of 
the relative incident kinetic energy. The rates of inelas- 
tic scattering into individual magnetic sublevels are sep- 
arated. The magnetic field is fixed at B = 100 G. As was 
expected from the analysis of adiabatic curves shown in 
Fig. the sublevel Aid = 3 is responsible for most of 
the inelastic losses. 
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FIG. 3: The elastic and inelastic collision rates are shown 
as functions of the relative incident kinetic energy. This plot 
separates the contributions to the total inelastic rate from 
different Zeeman levels. The contribution from the lowest 
Zeeman levels are not shown since they are very small. For 
this plot, we have included contributions from all higher L > 
partial waves in the entrance channel. The higher partial 
waves are important at energies above 100 ^iK. 

In view of experimental interest in other metastable 
alkaline-earth atoms, we have carried out a similar cal- 
culation for calcium atoms in the (4s4p)^P2 state and 
obtained a comparable result for the elastic and inelas- 
tic rates. The main difference is that the resonance for 
Ca2 is shifted toward larger magnetic fields. To as- 
sess whether any other species is likely to produce a 
more favorable elastic/inelastic ratio, we have completed 
a model calculation in which we artificially change the 
quadrupole moment, the nuclear mass and the disper- 
sion coefiicients (see Refs. 0, ^| for definitions of the 



quantities) . Although such a variation of the atomic pa- 
rameters modifies the results, the main result remains 
the same, namely an unfavorable elastic/inelastic ratio 
at submillikelvin coUisional energies. We point out that 
a similar result was recently obtained in the analysis of 
the collisions between two polar molecules jlQ] in an ex- 
ternal electric field. Thus, large inelastic rates apparently 
rule out the possibility of cooling atoms or molecules with 
large anisotropic multipole interactions down to submi- 
crokelvin temperatures using magnetic traps. 

A pos sible solution to this problem is to use an optical 
trap [I3 instead of a magnetic one. In this context, the 
situation is similar to the case studied in atomic cesium. 
Early experiments attempting Bose-Einstein condensate 
of Cs (see, for example [3) using magnetic trapping ap- 
proaches had shown that the degenerate regime can not 
be reached due to a rapid loss of atoms at short internu- 
clear distances. In the case of Cs, the losses were caused 
by a large spin-flip cross-section. Recently, the group of 
Grimm 17] utilized a purely optical trap of Cs atoms and 
successfully achieved a Cs Bose-Einstein condensate. 

Conclusion. We have calculated the elastic and in- 
elastic two-body collisional rates of metastable strontium 
atoms in the [(5s5p)^P2] state and metastable calcium 
atoms in the [{AsAp)^P2\ state. At low temperatures, 
< \^K, collisions in all low-field seeking states are domi- 
nated by inelastic processes. This appears to rule out the 
possiblity of achieving a Bose-Einstein condensate with 
these metastable atoms, using magnetic traps. 
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